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SUMMARY 


Results are presented of a theoretical investigation of an auto- 
matic control system with primary sensitivity to normal accelerations 
as used to control a specific supersonic canard missile. The acceler- 
ation control system consists of an integrating servomotor which 
receives its actuating signal from an accelerometer sensitive to the 
normal accelerations of the airframe being controlled. The servomotor 
operates the airframe control surfaces to obtain or maintain desired 
normal accelerations. 

The analysis is based largely on comparisons of normal acceleration 
transient responses obtained for various conditions of Mach number, 
altitude, static margin, and rate-of -pitch feedback. The results indi- 
cate that the use of rate-of-pitch feedback and a high static margin 
with accompanying increase in integrating- servomotor gain and rapid 
control-surface deflection results in a more rapid transient response 
and a lower steady-state attitude error due to aerodynamic out-of-trim 
moment . 

The acceleration control system appears to be a satisfactory method 
for obtaining longitudinal control of the supersonic airframe under con- 
sideration. This system has no directional space reference of its own, 
however, and its primary usefulness is therefore believed to be in con- 
junction with a homing seeker or with a guidance system which will pro- 
vide a directional space reference. 
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INTRODUCTION 


As part of the general research program for investigating various 
means of automatic stabilization, the Pilotless Aircraft Research 
Division of the Langley Aeronautical Laboratory has been conducting a 
theoretical analysis to determine the possibilities of using an auto- 
pilot primarily sensitive to linear accelerations for longitudinal 
stabilization and control of a supersonic canard airframe. Physically 
this type of control combines the use of a linear accelerometer and 
servomotor to obtain desired normal accelerations of the missile. Since 
the acceleration control system has no directional space reference of 
its own, its primary usefulness is believed to be in conjunction with a 
homing seeker. An autopilot of the type investigated in this analysis 
is small, lightweight, and simple to fabricate as compared to a 
displacement-plus-rate type of autopilot. Longitudinal control through 
the use of an accelerometer also has the advantage of eliminating the 
problem of free gyroscope drift. 

The analysis has been made for a specific supersonic canard missile 
configuration and is based mainly on flight conditions and stability 
characteristics anticipated as a result of previous flight tests of 
geometrically similar models. The results presented show the effects 
of the following conditions on the over-all performance of the autopilot- 
model combination: 

(1) Static margin and Mach number variation 

(2) The addition of a rate-of-pitch feedback control 

(3) Altitude variation 

(k) Aerodynamic out-of-trim moment 

(5) Variation of servo-gain constants 

(6) Mach number variation with fixed servo-gain constants 

(7) The use of an accelerometer displaced ahead of model center of 

gravity to generate both normal-acceleration error and rate- 
of-pitch feedback signal. 
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SYMBOLS 

normal acceleration of airframe in g units 
desired normal acceleration of airframe in g units 
voltage proportional to n G , volts 

voltage proportional to , volts 

voltage error, volts = V n ^ - V nQ ^ , or normal-acceleration 
error, g units (e = n^ - 

acceleration due to gravity, 32.2 feet per second per second 

pitch angle measured from the horizontal, degrees 

first derivative of pitch angle with respect to time, 
degrees per second (d0 o /dt) 

angle of attack, degrees 

first derivative of angle of attack with respect to time, 
degrees per second (dxx/dt) 

flight-path angle, degrees (7 = Q - a) 

first derivative of 7 with respect to time, degrees per 
second (d7/dt) 

canard control-surface deflection, degrees ^5 = 6g - 6^ y 

control- surface deflection due to integrating servo, degrees 

control- surface deflection due to rate servo, degrees 

value of control-surface deflection which counterbalances 
out-of-trim moment, degrees 

integrating servo gain constant, radians per second per g 
rate-servo gain constant, radians per radian per second 
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steady-state proportionality constant between voltage and 


time, seconds 
Mach number 

stability axis which passes through center of gravity and 
is perpendicular to vertical plane of symmetry 

moment of inertia about Y-axis, slug-feet square 

mean aerodynamic chord, feet 

model center of pressure 

static margin, fraction of c 

velocity, feet per second 

mass, slugs 

• _ 

0 C 

dynamic pressure, pounds per square foot; or — (when 
used as a subscript) 

wing area, square feet 


acceleration, volts per g 



lift coefficient 



pitching-moment coefficient 
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c mq /% if 

i \[^L 

oo angular frequency, radians per second 

D differential operator (d/dt) 

s Laplace transform variable corresponding to differential 

operator 

KG system or component transfer function; may be expressed as 

a function of jo>, D, or s 

AR magnitude of KG(joo) 

PA phase angle of KG(jco) 

R Routh's discriminant 


A,B,C,D,E,F 


coefficients of the quintic characteristic equation 



DESCRIPTION OF THE PROPOSED ACCELERATION CONTROL SYSTEM 


The block diagram of the proposed acceleration control system with 
the airframe compensated by rate-of-pitch feedback is shown in figure 1. 
The voltage input v ni of the system is made proportional to a desired 

normal acceleration n.j_ of the airframe. The accelerometer produces a 
voltage V no proportional to the existing normal acceleration n Q of 
the airframe. The proportionality constant between desired acceler- 

ation n^ and input voltage V n -^ is the same as the proportionality 
constant between existing acceleration no and acceleration volt- 
age Vtiq- Then if existing acceleration is not equal to desired acceler- 
ation an error signal € excites the integrating servo. The integrating 
servo produces a control -surface deflection 5 g which is proportional 
to the integral of the error signal. This control surface deflection 
causes the airframe to turn in the proper direction to produce a normal- 
acceleration signal which tends to cancel the error signal. In the 
steady-state condition the error signal is zero, but the integral of the 
error, and hence the control-surface deflection, is not necessarily zero. 
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The effect of the rate servo can he thought of as a modification of 
the airframe response since its primary effect is to increase the 
damping ratio of the airframe. 

Proposed configurations for the various elements, other than the 
airframe, are presented in the following paragraphs. 

Rate servo .- The proposed rate servo consists of a rate gyro and 
hydraulic servo combined as shown schematically in figure 2(a). Fig- 
ure 2 illustrates one possible mechanical arrangement, many variations of 
which are possible. If time lags and the effect of inertias are ignored, 
the control deflection 5 -r is proportional to a rate input 0. For the 
present analysis, the rate servo was considered to be a s ingle- degree- 
of-freedom, second-order system having a natural frequency of 88 radians 
per second and a damping ratio of 0.5- This is believed to be a con- 
servative representation of the dynamic effects of an actual rate gyro 
plus servomotor. 

Integrating servo .- A proposed integrating servo is shown schemati- 
cally in figure 2(b). Since there are moving masses in this servo the 
response cannot be instantaneous as implied by the transfer func- 

tion, — = _1. However, the assumption made herein i6 that the fluid- 
supply pressure is high with no limit on the rate of flow, and the 
masses of the moving parts (including moving fluid) are small. Since 
such a servo is fast acting, the dynamic effects can be ignored for 
operation at low frequencies. 

Accelerometer dynamics .- The term "accelerometer dynamics" as used 
here includes any dynamic effect between the normal acceleration n Q 
of the airframe and the electrical accelerometer signal V no . Most 
accelerometers having a range suitable for the present application have 
a natural frequency which is sufficiently high so that the dynamic 
effects of the accelerometer may be neglected; however, because of noise 
signals picked up by the accelerometer, the accelerometer signal may 
have to be filtered electrically. The dynamic effects of such a filter 
are included in the accelerometer dynamics discussed here. 

For the major part of the present investigation, the effect of 
accelerometer dynamics was neglected; however, its effect was investi- 
gated for one set of conditions, and the results obtained are shown 
later. For this set of conditions, the following transfer function was 
used to represent the accelerometer dynamics. 

Vn o = 246 50 

n o D 2 + 3l4D + 24650 
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This transfer function represents a second-order system having a natural 
frequency of 25 cycles per second and damping ratio of 1.0. At 5 cycles 
per second the amplitude ratio is 0. S6 and the phase lag is 23°, and 
at 80 cycles per second the amplitude ratio is 0. 1. This transfer func- 
tion was used because it is believed that these phase and attenuation 
characteristics are representative of those to be expected from an actual 
accelerometer and filter suitable for the present application. 


METHOD OF ANALYSIS 


The analysis is mainly concerned with obtaining the transient accel- 
eration response (no] to a unit-step acceleration input (ni] . These 
responses are then used to determine the effect of variation of the 
automatic -control system gain constants and the airframe aerodynamic 
parameters. To facilitate analysis, the airframe and autopilot com- 
ponents can be represented by transfer functions which can be combined 
in block diagram form. 


Component Transfer Functions 

Servomotor block . - This block represents an integrating servo. If 
the dynamics of the servomotor and control surface are neglected, this 
block can be drawn as : 


§ 

K 1 

6 S 


K l °i * IT 



Airframe block. - The forms of the equations of motion for constant 
speed and disturbances from level flight are: 

(jjf 1)2 ' c ”q #? D ) e ° " (°"ta + m D ) ° “ C H6 5 

(f E )'o - (S D + %)“ * 

where the stability derivatives are expressed in radian measure. 


8 


NACA RM L51D23 


Q 

Solution of these equations for gives the airframe transfer 

function, which can be represented in block diagram form as: 


O ^ 

ir r a(D + b) _ 

e 0 

^ 

K 2 G 2 D(D 2 + e D + f) 



where 


a = 


Cm s(cF 


C L&( Cn \x it 



C m Pl C I w “ C m n C LR 

°i„ (% + c ^)l? 


° I tt( Cm CL 2v) + Cm a(qsj 

fi3L\ « 

\qScj \qS/ 



The np/Rp block .- The transfer function n O/ /0 O is obtained from 
the relation 

Uq _ n 0 /S 

®o/^ 

where 

B2 _ v Z 

5 32.2 5 


7 

and the function — 

o 

equations of motion. 


is obtained by substituting a + 7 
The solution for g- is then 


7 

5 
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^o 


in the 
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For the case where 
units per radian can be 


C-Lg = 0^ the transfer function n o /0 o 
represented in block-diagram form as: 


in g 



KiG? - ^ 

no 

^ 

3 3 a (D + b) 



where 

CL. C T 

h = I m s hi 

32.2 / 1 y \ / mV\ 

U Sc/ \qS] 

and a and b are previously defined. 


Rate-control block . - The rate control can be represented in block- 
diagram form as: 


k r°r 


jnhh k r d 
D 2 + 88D + 77kh 


Jk. 


This transfer function has been used previously in conjunction 
with a displacement autopilot in reference 1. The static sensitivity 
( k r) has the physical significance of being the magnitude of steady- 
state control-surface displacement 5p, resulting from a unit rate 
input 9. 


Accelerometer block .- For the major part of this analysis, the 
dynamics of the accelerometer in^the outside loop have been neglected. 

That is, the transfer function is assumed to be equal to unity. 

In determining the effect of including the accelerometer dynamics, it 
was assumed that the accelerometer and filter could be represented by 
the following block: 



The gain was taken as unity since variations of have the same 
effect on the responses to an n-^ input as do variations of servo gain 
Kq, as can be derived by using the block diagram shown in figure 1. 
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Method of Obtaining Transient Responses 


The system analyses presented herein are based mainly on the tran- 
sient responses of the system to a unit step input. Although such dis- 
continuities in the input may never occur in practice, the responses to 
a unit step input are of value in making a comparison of a system’s 
performance under various conditions. 

The majority of the transient responses presented herein are square- 
wave responses obtained from plots of the closed-loop frequency response 

— (jo>) by the use of the following series: 
n i 


°o(t> = (f- 


+ £ 


<n=0 


z 


m=l,3,5. 


AR 

m 


sin 


fmuqt + PA) 


where n Q (t) is the response to a square-wave input. If the funda- 
mental frequency of the square wave is made low enough so that the 

system transients die out in one-half period of the square wave, the 
square-wave reponse is essentially the response to a series of steps. 
Twelve terms of the above series were used and summed by a Fourier 
synthesizer. This instrument was developed by the Massachusetts 
Institute of Technology, and the mechanical and technical aspects of a 
similar instrument along with a description of the equipment and a 
derivation of the preceding series can be found in reference 2. 

A graphical method employing the general theories of servomechanism 
analysis, as outlined in reference 3, was used to obtain the closed-loop 
frequency response prior to obtaining each transient. With the assump- 
tion made previously = 1^ , the block diagram for the combination of 

automatic control system and model reduces to: 
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This analysis is mainly concerned with the transfer function 
from which the transient response n Q (t) to a unit step n^ can be 
obtained. As mentioned previously, a graphical method as described in 
reference 3 (chapter 8) is employed in obtaining Since the block 

diagram contains two feedback loops, the solution must be handled in 
two steps. The first step consists of obtaining an Lm-Angle plot of 

the open-loop response -£■, on which the M-N contours are super- 
imposed. Then a satisfactory (but not necessarily optimum) adjustment 
of the rate servo gain constant Kp can be made and the closed-loop 

8p 

response — obtained. The second step combines the inner loop with 

the remaining components of the block diagram to obtain the Lm-Angle 
plot of the open-loop response — on which the M-N contours are also 
superimposed. The integrating servo gain constant is adjusted on 

this plot. The significance of this adjustment is that it fixes the 
peak amplitude ratio of the closed- loop ^(jco). In most cases, the 
value of K x was chosen, as suggested in references 3 and k, so that 

>o>) = 1.3 

n i max 

RESULTS AND DISCUSSION 


The results which follow show the effect of varying the airframe 
static margin, flight conditions, and autopilot gain constants on the 
over-all performance of the combination of the automatic-control system 
and model using a specific supersonic airframe. A photograph and plan- 
view sketch of the airframe are shown in figure 3» Flight tests of 
geometrically similar airframes have previously been conducted, the 
results of which are presented in references 6, and 7- The estimated 
and measured longitudinal derivatives given in table I are based on 
reference 6, in which the model center of gravity was 73- 53 inches behind 
the nose. The measured derivatives for intermediate static-margin values 
were taken directly from this reference, and the values given for small 
and large static margins were obtained from the estimated changes of the 
measured derivatives due to shifts in the airframe center of gravity. 

The static margin varies with Mach number in each static-margin category 
listed in table I. The values given for intermediate static margin, for 
example, vary from 0.333c to 0.29kc as M varies between 1.0 and 2.0. 
This variation is due to shifts in the airframe center of pressure wi+n 
Mach number. The variations of the standard atmosphere flight condit ons 
used in this analysis are given in table II. 
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Effect of static margin and Mach number variation .- The transient 
responses (n Q to a unit step acceleration input) shown in figures 4, 5> 
and 6 are for small, intermediate, and large static margins, respec- 
tively. On these figures, the value of the rate-servo gain constant 
(Kr) was chosen such that the value of the peak amplitude ratio of the 

transfer function ^rr(jo)) was equal to 1.3 and the value of the 


4 (Ja) 

integrating-servo gain constant 


was chosen such that the peak ampli- 


tude ratio of the closed-loop frequency response ^(joi) was also equal 

to 1. 3. An examination of figures 4, 5, and 6 reveals that the response 
time (the time required for the output transient to reach and remain 
within a given percent of steady state) is decreased by increasing either 
static margin or Mach number. The most rapid responses obtained are 
shown in figure 6 which is based on sea- level flight of the model with 
large static margin. 


Effect of rate-of -pitch feedback . - Some of the responses shown in 
figures 4, 5> and 6 are for Kp = Oj setting Kp equal to zero has the 
effect of removing the rate-servo block (see fig. l) . An example of the 
effect of rate-of-pitch feedback is best shown by referring to figure 5 
where the system transient response for Kp =0 is shown along with 
responses including 9 feedback for each Mach number. An examination 
of the responses shown in this figure indicates that including 6 feed- 
back in this type of control system has the effect of increasing the 
damping of the missile and allows an increase in the integrating-servo 
gain constant K^, which results in a faster response. 

As mentioned in the section METHOD OF ANALYSIS, the transient 
responses were obtained by the use of a Fourier synthesizer. A com- 
parison between a Fourier synthesizer result and a transient response 

calculated by the methods of Laplace (references 3 and 8) is made in 

figure 7 for the case of sea- level flight at M = 1.6 and with 
SM = 0.294c. Other comparisons of the results obtained by these two 
methods have been made and agreement between the two methods is con- 
sidered sufficient to justify use of the Fourier synthesizer to obtain 

the transient responses. 


Effect of altitude variation . - Except for altitude and servo- gain 
adjustments, the responses shown in figure 5(c) are for the same con- 
ditions as those presented in figure 8 where the n Q transient responses 
to a unit step acceleration input are based on flight at 10,000 feet 
and 40,000 feet. The values of K-l and Kp used in figure 8 (except 

for the dotted curve of figure 8(a)) were chosen such that the peak 

S-r n Q , \ 

amplitude ratios of the transfer functions (j«>) and ^ ^ were 
equal to 1.3. As is shown, flight at 10,000 feet produces a somewhat 
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slower response than that obtained at sea level while flight at 
40,000 feet shows a considerable increase in the normal acceleration 
response time. A comparison of figure 5(c) and figure 8(a) indicates 
that the value of these servo gain constants did not change appreciably 
between sea level and 10,000 feet. The dotted curve of figure 8(a) is 
the response at 10,000 feet with rate-of -pitch feedback using the same 
servo gain adjustments that were used to obtain the sea-level response 
of figure 5(c). Since these responses do not differ greatly, it is 
believed that a fair approximation of the behavior of the airframe 
between sea level and 10,000 feet is obtained by basing the analysis 
entirely on sea- level flight. 

Effect of accelerometer dynamics .- As mentioned previously in the 
description of the proposed acceleration control system, the acceler- 
ometer dynamics can be represented by the transfer function 



2^650 

D 2 + 3 i 4 d + 2^650 


In figure 1, this transfer function is labeled "accelerometer dynamics" 
in the outer feedback loop. In figure 9 , a comparison is made of the 
response obtained by using the proposed acceleration control system with 
and without the accelerometer dynamics included in the outer feedback 
loop. The curves shown in figure 9 are based on sea-level flight at 
M = 1.6 and with SM = 0.29^c, and the results indicate that a slight 
increase in response time and period of the transient oscillations is 
obtained with the inclusion of accelerometer dynamics. On this basis, 

Vn^ 

using the simplifying assumption that = 1 for the major part of 

this analysis seems to be justified. 

Effect of aerodynamic out-of-trim moment . - In the automatic control 
system being investigated, an aerodynamic out-of-trim moment can be 
represented by an equivalent control- surface deflection, to be denoted 


Ik 
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"by 5 o1 _ . This control- surface deflection can be represented in the 
system block diagram as an additional input as follows: 


6ot 



The transfer function ^2— can be derived by setting n* = 0 and 

°ot 

employing the relations 

n o = - c 

and 

5 = 5 S - 5 R + 5 ot 

Using the general theories of servomechanism analysis, the transfer 

function is: 

5 ot 

E°_ = f(s) (X) 

§ot 1 + t KgGgK-pG-j^ 

The steady-state error is then obtained by applying the final-value 

theorem (see reference 3> chapter 3) "to equation (l) as follows: 


lim n Q (t) = 
t 



00 


lim s 
s — 0 
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from which it is found that (steady state) is always zero for the 

Oot 

control system under consideration. 

e r 


Similarly, the transfer function 


'_ o 

Sot 


can be derived as 


^£2 

&ot 1 + K 1 G 1 K 2 G 2 K3G3 + k 2 g 2 k r g ]r 


( 2 ) 


Then by applying the final-value theorem, —— (steady state) is found to 

Oot 

be finite, the value of which is a measure of the effect of out-of-trim 
moment . 


The results of the out-of-trim investigation conducted herein are 

presented in table III where the values of (steady state) are tabu- 

hot 

lated for the automatic stabilization system with and without rate-of- 
pitch feedback. The results show that the system with rate feedback 
produces less 0 o error due to 5 0 -^ for any set of comparable con- 
ditions. It is also shown that the magnitude of the error increases 
with increasing Mach number and decreases with increasing static margin. 
These results indicate that the use of rate feedback and high static 
margin will keep the Q q error due to out-of-trim moment at a minimum. 

To summarize, the principal effect of an aerodynamic out-of-trim moment 
is to cause a steady-state error in pitch angle 0 O , with no steady-state 
error in normal acceleration. 

Application of Routh's discriminant to the stability analysis . - 
An application of Routh’s criterion for the case of sea-level flight at 
M = 1.6 with SM = 0.29^c is presented in figure 10. The characteristic 
equation of the closed-loop transfer function ^j(s) with rate-of-pitch 

feedback is of fifth degree. The conditions for complete stability for 
a system having a characteristic equation of fifth degree are derived in 
reference 9, and applications of Routh’s criterion including the neces- 
sary and sufficient conditions for complete stability are given in 
reference 10. 

In conducting the analysis based on Routh's criterion as presented 
herein, the coefficients of the characteristic stability equation were 
expressed as functions of K-^ and K R . The curve presented in fig- 
ure 10 was then obtained from the condition for neutral oscillatory 
stability 


R = (BC - AD) (DE - CF) - (BE - AF) 2 = 0 
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Any combination of the values of Kp and Kp falling above the neutral 
oscillatory stability boundary in figure 10 produces oscillatory insta- 
bility. The lower stability limit of the integrating servo gain Kp 
was determined from the coefficient F which is a constant multiplied 
by Kp. Therefore, since an unstable root would exist if Kp assumed 
a negative value, the lower limit for Kp is zero. Since the other 
conditions for complete stability as given in reference 10 fall outside 
but indicate stability in the direction of the stable region of fig- 
ure 10, any set of values of Kp and Kp falling within this region 
will produce a stable system. 

Effect of varying rate-servo gain constant (Kp).- A locus of points 
based on six values of Kp and for which Kp was adjusted to make the 
peak- amplitude ratio of the transfer function j^-(jao) equal to 1.3 is 

shown in figure 10. Closed-loop frequency responses ^-(jo)) based on 
values of Kp and Kp which correspond to the points of this locus and 
based on the same flight conditions used in figure 10 are presented in 
figure 11. In figure 12, the transient responses (n Q to a unit step np 
obtained from each frequency response are shown. The results given in 
figure 11 indicate that the resonant frequency peak occurs at a higher 
frequency as Kp is decreased, and as the value of Kp approaches zero 
a dip, or bucket, appears in the lower frequency range of the amplitude- 
ratio curves. In figure 12, a slight decrease in the transient response 
time is shown as Kp decreases; however, for the case of Kp = 0, a 
slowly rising oscillatory transient response, resulting principally 
from the low frequency characteristics of the frequency response, is 
obtained as shown by the final plot. The foregoing data indicate that 
values of Kp in the range 0.062 to 0.049 yield satisfactory responses 
in that they are rapid and well damped, although the exact adjustment 
of Kp does not seem to be critical if Kp is adjusted to obtain a 
slight overshoot of the n Q transient response. 

Effect of varying Mach number with servo gains fixed. - The effect 
of fixing the values of the servo gain constants and varying the Mach 
number is shown for intermediate and large static margins in figures 13 
and 14, respectively. The curves presented in these figures are based 
on sea-level flight at four different Mach numbers and show the normal 
acceleration responses obtained when Kp and Kp are fixed at the 
values shown in figure 15 at M = 1.6. 

The variations of Kp and Kp with Mach number shown in figure 15 
are based on the values of these gain constants used previously in 
obtaining the responses shown in figures 5 and 6 where the peak ampli- 
tude ratio of ^(jo>) was set at 1.3- It is shown in figures 13 and l4 
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that , when is adjusted to give 


&<*»» 


max 


= 1. 3 at M = 1.6, 


values of M less than 1.6 give more stable responses and the n Q 
transient response becomes more oscillatory as M is increased above 1.6. 
The oscillatory response obtained for M = 2.0 in figure l4 indicates 
that the effect of varying Mach number with servo gains fixed is greater 
with large static margin. 


The use of an auxiliary control to vary the servomotor gains with 
Mach number to conform with the values given in figure 15 would reduce 
the effect of Mach number variation. Since it has already been shown 
that the exact adjustment of K R is not critical if K-^ is adjusted to 
obtain a slight overshoot of the n 0 transient response, the design of 
a gain- varying device could be simplified by basing it only on the vari- 
ation of K]l with Mach number. 

Transient responses of 0 O , a, 7, and 6 to a unit step acceler- 
ation input .- The results presented in figure l6 are based on sea- level 
flight at M = 1.6 for two values of static margin. Figure 16(a) is 
for SM = 0.29^5, whereas figure 16(b) is for SM = 0.564c. Normal- 
acceleration transient responses to a unit step acceleration input for 
the same conditions as used in obtaining the responses of figure 16 have 
previously been shown in figures 5(c) and 6(c). It can be seen in fig- 
ure 16 that the angle of attack approaches a steady-state value of 
approximately 0.23° per g normal acceleration for either value of static 
margin, and after 0.3 second the steady-state rate of change of pitch 
angle and flight-path angle is approximately 1° per second per g. An 
examination of the control-surface-deflection (&) responses presented 
in figure 16 shows that , as would be expected, approximately twice as 
much steady-state 6 is required per g normal acceleration with the 
larger static margin; namely, 0.49° per g for SM = 0.564c and 
0.255° per g for SM = 0.294c. 

A further examination of figures 16(a) and l6(b) shows that more 
rapid missile responses are obtained with high static margin. This 
result is obtained because the increased aerodynamic stability of the 
airframe allows a higher integrating- servo gain K^, which produces a 
more rapid change in control-surface deflection. 


Gravitational influence. - The results of this analysis can be 
applied either to pitch or to yaw due to the symmetry of the airframe 
in these planes. However, in pitch, because of the 1 g acceleration 
due to gravity, there exists a gravitational influence on the acceler- 
ometer which varies as cos Q 0 . This influence has been neglected 
throughout this analysis. Its omission is not believed to affect the 
transient responses seriously since the total variation of Q q during 
the transient response time is small, as can be seen in figure 16. The 
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effect of this gravitational influence on the trim condition is not 
serious because in practice the trim condition will be continuously dic- 
tated by a homing seeker or other guidance system. 

System responses obtained with accelerometer placed ahead of center 
of gravity to generate rate-of -pitch feedback signal .- In practice, it 
may be desirable to include only one servo in the system and add the 
normal-acceleration and rate-of-pitch signals electrically or by some 
other method. Mr. H. D. Garner of the Instrument Research Division at 
the Langley Laboratory has suggested that the rate-of-pitch feedback 
may be obtained by mounting the normal accelerometer ahead of the center 
of gravity of the airframe so that the accelerometer will be sensitive 
to angular acceleration §’ as well as normal acceleration iiq. Since 
the accelerometer signal is fed through the integrating servo, the com- 
ponent due to angular acceleration is effectively an angular-rate feed- 
back or rate-of-pitch feedback. 

The rate-of-pitch feedback gain is then determined by the distance 
from the center of gravity of the airframe to the normal accelerometer. 
With static margin of 0.294c and M = 1.6, to obtain the response of 
figure 5(c) (with Kp = 0.11 and K R = O.062), the normal accelerometer 
must be mounted 18.1 feet ahead of the airframe center of gravity. With 
static margin of 0.564c and M = 1.6, to obtain the response of fig- 
ure 6(c) (with Kp = 0.21 and K R = 0.047), the normal accelerometer 

must be mounted 7.2 feet ahead of the airframe center of gravity. The 
size of the airframe used in this investigation limits the distance 
between airframe center of gravity and normal accelerometer to a maxi- 
mum of about 5 feet. 

A theoretical investigation of the system response with the normal 
accelerometer mounted 5 feet ahead of the airframe center of gravity 
with static margin of 0.564c and Kp = 0.09 (effective K R = 0.0l4) was 
made for Mach numbers of 1.0, 1.2, 1.6, and 2.0. The results are shown 
in figure 17. It can be seen from figure 17 that this location of the 
accelerometer does not provide enough effective pitch-rate feedback to 
damp the system satisfactorily. 

To produce sufficient damping by this method, it would be necessary 
to extend the nose of the model 2 feet or more. An alternative is to 
use two normal accelerometers, one mounted a distance Z ahead of the 
airframe center of gravity and one mounted on the airframe center of 
gravity. The signals from these accelerometers are subtracted before 
being fed to the integrating servo. The pitch-rate feedback gain and 
the acceleration feedback gain can then each be adjusted independently 
for any distance Z by adjusting independently the two accelerometer 
gains and the integrating servo gain. It is believed that no difficulty 
will be encountered due to slight mismatching of the dynamics of the 
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two accelerometers since the accelerometer natural frequencies will be 
considerably higher than the resonant frequency of the system and since 
the high frequency components of the accelerometer signals will be 
attenuated by the integrating servo. 

Comparison of acceleration control system with an attitude control 
system . - A reasonable way to compare two control systems to be used for 
missile guidance is to compare their effectiveness in obtaining rapid 
changes in flight-path direction without producing excessive normal 
accelerations of the airframe. This comparison has been made between 
the acceleration control system analyzed herein and the attitude control 
system analyzed in reference 1; the results axe shown in figure 18. The 
results for the acceleration control system and for the attitude control 
system were compared at a Mach number of 1.6 and 1.8, respectively, 
because these results were available. The static margins used were con- 
sidered to produce the most satisfactory results for each of the two 
systems, that is, intermediate static margin for the acceleration control 
and small static margin for the attitude control. In figure l8, the com- 
mand input 6± to the attitude control system is a one-degree step. 

As is shown, this produces a peak normal acceleration of 4.1g and results 
in a steady-state change in flight-path angle y of 1°. The y tran- 
sient has died out to within 5 percent of its final value in approximately 
0.86 seconds. The command input to the acceleration control system is a 
square pulse. The magnitude of the square pulse was chosen so as to pro- 
duce a peak normal acceleration of 4. lg and the time duration was chosen 
to produce a steady-state change in flight-path angle y of 1°. As 
shown in figure 18 , the y transient has died out to within 5 percent 
of its final value In approximately 0.45 second or approximately one- 
half the time required by the attitude control system. This result 
indicates that changes in flight-path direction may be obtained more 
rapidly with the acceleration control system. 

The areas under each of the normal-accelerat ion-response curves 
shown in figure 18 are approximately equal. Since the normal acceler- 
ation is proportional to angle of attack, there is no apparent increase 
in velocity loss due to drag in obtaining the shorter y response time 
with the acceleration control system. The comparison made here is based 
on the available data for the attitude control system. It is believed 
that more cases and other variables would have to be considered before 
a general conclusion as to the relative merit of either system can be 
made. 
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CONCLUSIONS 


The acceleration control system is primarily sensitive to normal 
accelerations of the airframe and actuates the control surfaces through 
the use of an integrating servomotor to obtain desired normal acceler- 
ations. An automatic control system with primary sensitivity to linear 
acceleration appears to he a satisfactory method for obtaining longi- 
tudinal control of the supersonic airframe under consideration. The 
acceleration control system, however, has no directional space reference 
of its own; therefore it is believed that the primary usefulness of such 
a system is in conjunction with a homing seeker or with a guidance system 
which will provide a directional space reference. The conclusions reached 
as a result of the analyses presented herein, based on a specific super- 
sonic airframe configuration, are as follows: 

1. On the basis of a comparison of the normal-acceleration tran- 
sient responses presented herein for various Mach numbers, static margins, 
and altitudes, it can be concluded that: 

(a) Increasing the airframe static margin produces more rapid 
transient responses, which may necessitate the use of a fast-acting 
servomotor. 

(b) Including rate-of-pitch feedback in the automatic control 
system has the effect of increasing the damping of the missile and 
allows an increase in the integrating- servo gain constant resulting 
in a more rapid response. The exact adjustment of the rate-of-pitch 
feedback control gain constant does not seem to be critical if the 
integrating-servo gain is adjusted to obtain a slight overshoot of 
the normal-acceleration transient response. 

(c) More rapid transient responses can be obtained with higher 
Mach number. 

(d) Flight at altitude produces slower responses than those 
obtained at sea level; however, a fair approximation of the behavior 
of the system between sea level and 10,000 feet is obtained by 
basing the analysis entirely on sea-level flight. 

2. When the integrating-servomotor gain constant is adjusted so 
that the peak amplitude ratio of the system-closed-loop transfer func- 
tion is 1. 3 at a Mach number of 1.6, flight at lower Mach numbers yields 
more stable transient responses and the system tiansient response 
becomes more oscillatory as the Mach number is increased above 1.6. The 
effect of varying Mach number with fixed servomotor gain constants is 
greater with large static margin. 
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3. The principal effect of an aerodynamic out-of-trim moment is to 
cause a steady-state error in pitch angle with no steady-state error in 
normal acceleration. The use of a rate-of -pitch feedback control and 
high static margin tends to minimize the steady-state error in pitch 
angle due to an aerodynamic out-of-trim moment of the airframe. 

4. A comparison between the acceleration control system and an 
attitude control system indicates that changes in flight-path direction 
may be obtained more rapidly with the acceleration control system with 
no apparent increase in velocity loss due to drag. 

5. Theoretically, effective rate-of -pitch feedback may be obtained 
in a system employing an integrating servomotor by using a normal 
accelerometer mounted ahead of the model center of gravity; thus the 
necessity of a rate gyro in the control system may be eliminated. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 

LONGITUDINAL DERIVATIVES 


(All derivatives in radian measure; Iy = 31*3 
m = 5-05 slugs; C = 1.776 ft; S = 4. 1 


slug-ft^; 
sq ft] 


Mach 

number 

SM 

(fraction of c) 

Cinq 

Cm a 

c m§ 

Cn >d 

Cl 8 

°I<X 


Estimated derivatives for small static margin 
(static margin = 2 inches at M = 1.6) 

1.0 

0.126 

-6.55 

-0.406 

0.797 

-0.730 

0 

3.22 

1.2 

.134 

-6.96 

-.404 

.820 

-770 

0 

3.02 

1.6 

.094 

-6. 68 

-.245 

.702 

-740 

0 

2.61 

2.0 

.042 

-6.61 

-.104 

.573 

-.730 

0 

2.46 

Measured derivatives for intermediate static margin 
(static margin = 6.27 inches at M = 1.6) 

1.0 

0.333 

-6.97 

-1.07 

0.797 

-0. 774 

0 

3.22 

1.2 

• 339 

-7.48 

-1.025 

.820 

-.831 

0 

3-02 

1.6 

.294 

-7.22 

-.763 

.702 

-.802 

0 

2.61 

2.0 

.249 

-6.39 

-.613 

• 573 

-.717 

0 

2.46 


Estimated derivatives for large static margin 
(static margin = 12 inches at M = 1.6) 

1.0 

0.601 

-8.13 

-1.94 

0.797 

-0.903 

0 

3.22 

1.2 

. 606 

-8.88 

-1.83 

.820 

-.987 

0 

3.02 

1.6 

.564 

-8.47 

-1.47 

.702 

-.941 

0 

2.61 

2.0 

.516 

-7-43 

-1.27 

• 573 

-.826 

0 

2.46 
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TABLE II 

VARIATION OF FLIGHT CONDITIONS 


Altitude 

(ft) 

Mach 

number 

q. 

(lb/ sq ft) 

V 

(ft/sec) 

Sea level 

1.0 

l48l 

1116 

Sea level 

1.2 

2132 

1339 

Sea level 

1.6 

3791 

1785 

Sea level 

2.0 

5980 

224-5 

10,000 

1.6 

2602 

1721 

40,000 

1.6 

702 

1553 
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TABLE III 

STEADY- STATE 0 Q ERROR DUE TO S ot 


Mach 

Number 

SM 

(fraction of c) 

Altitude 

(ft) 

6 

(steady state) 

°ot 

Without rate-of-pitch feedback 

1.0 

0.126 

Sea level 

1.62 

1.2 

• 134 

— do 

1.45 

1.6 

.094 

—do 

2.96 

2.0 

.042 

do 

5.86 

1.0 

• 333 

do 

0.8 

1.2 

• 339 

do 

.8 

1.6 

.294 

do 

• 9 

2.0 

.249 

do 

• 96 

1.6 

.564 

do 

• 53 

1.6 

.294 

10,000 

• 93 

l. 6 

.294 

40,000 

.69 


With rate-i 

of -pitch feedback 

1.6 

0.094 

Sea level 

1.08 

2.0 

.042 

—do 

2.15 

1.0 

• 333 

—do 

0.13 

1.2 

• 339 

—do 

.16 

1.6 

.294 

—do 

• 17 

2.0 

.249 

—do 

.21 

1.0 

.601 

do 

• 073 

1.2 

. 606 

—do 

• 075 

1.6 

.564 

do 

.085 

2.0 

. 516 

do 

.096 

1.6 

.294 

10,000 

.16 

1.6 

.294 

40,000 

.20 
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Figure 1.- Block diagram of the proposed acceleration control system. 
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Figure 2.- Schematic diagram of proposed rate servo and integrating servo 
illustrating a possible mechanical arrangement. 
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(a) Photograph of model configuration. 



Figure 3*“ Supersonic missile research model configuration. 
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(b) Plan-view sketch of model configuration. 


Figure 3.- Concluded. 
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Figure 4.- Normal acceleration (no) responses to a unit step acceleration 
input (ni) for various Mach numbers based on sea- level flight of the 
model with small static margin. Ki is the integrating- servo gain 
constant and Kr is the rate-servo gain constant. 
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(a) M=/.0 ; SM = 0.333c (b) M=/2 ; SM = 0.339c 




Cc) M=/6 ; SM=0.£94c (d)M=£.0 ; SM =b.£49c 


Figure 5.- Normal acceleration (no) responses to a unit step acceleration 
input (np) for various Mach numbers based on sea- level flight of the 
model with intermediate static margin. Kp is the integrating- servo 
gain constant and Kr is the rate-servo gain constant. 
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Figure 6.- Normal acceleration (no) responses to a unit step acceleration 
input (ni) for various Mach numbers based on sea- level flight of the 
model with large static margin. Ki is the integrating- servo gain 
constant and Kr is the rate-servo gain constant. 
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t, sec 


Figure 7>- Comparison of a transient response obtained using the Fourier 
synthesizer with one obtained using the methods of Laplace. Both 
curves are based on sea-level flight at M = 1.6 and SM = 0.294c. 

Kp = 0.11, K r = 0.062. 
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(b) 40,000 feet. 

Figure 8.- The effect of flight at altitude on the normal acceleration 
(n Q ) responses to a unit step acceleration input (ni) . Responses 
are based on flight at M = 1.6 and SM = 0.29^5. 
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Figure 9-- Comparison of the closed-loop transient responses for the 
system with and without the dynamics of an accelerometer included 
in the outer feedback loop. Both curves are based on sea-level 
flight at M = 1.6 and SM = 0.294c. (Kr = O. 062 ) 


I/jfegraP/rp serro go//?, /<, , radians per sec/g 
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Figure 10.- Application of Routh's discriminant to the characteristic 
equation of the closed-loop transfer function with the coefficients 
expressed as a function of Kp and Kr and based on sea-level 

flight at M = 1.6 and SM = 0.29^c. 
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Figure 11.- Closed- loop frequency responses 




n 


for various values 


of Kr and with values of K\ such that the peak amplitude ratio 


of 




M = 1.6 and 


1.3* All responses are based on sea-level 
SM = 0.29^c. 


flight at 


40 


NACA RM L5XD23 






Figure 12.- Normal acceleration (n 0 ) responses to a unit step acceleration 
input (ni) for various values of Kr and with values of Kp such that 


the peak amplitude ratio of 
sea-level flight at M = 1.6 



1.3. All responses are based on 


and SM = 0.29^-c. 
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(a) M=W j SM= 0.333 c 




Cc) M=/.6 ; SM= 0.294c 



Figure 13«- Normal acceleration (no) responses to a unit step acceleration 
input (ni) as obtained for fixed values of integrating- and rate- servo 
gain constants with Mach number as a variable. Responses are based on 
sea- level flight of the model with intermediate static margin. = O. 1-1 
and Kr = 0.062. 
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(d) M=2.0 ,SM= O.S/6 c 


Figure l4.- Normal acceleration (no) responses to a unit step acceleration 
input (np) as obtained for fixed values of integrating- and rate-servo 
gain constants with Mach number as a variable. Responses are based on 
sea- level flight of the model with large static margin. Kp = 0.21 

and Kp = O.OV 7 . 
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(a) SM = 0.29^-c; Kj = 0.11; Kp = 0.062. 

Figure 1 6 .- Transient responses of 0 Q y a, 7, and 5 to a unit step 
acceleration input based on sea- level flight at M = 1.6. 
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(b) SM = 0.564c; Kj_ = 0.21; Kr = 0.047. 
Figure 1 6 .- Concluded. 
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(a) Closed-loop frequency responses ^p(jcn). 


Figure 17*- System frequency and transient responses obtained for various 
Mach numbers with an accelerometer placed five feet ahead of the model 
center of gravity to generate both normal-acceleration error and rate- 
of-pitch feedback signal. Kp = 0.09; effective Kr = 0.0l4. 
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(b) Normal-acceleration (n 0 ) transient responses to a unit step acceleration 

input (n-jJ . 


Figure 17 •- Concluded. 
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Figure l 8 .- Comparison of 7 response times for acceleration control system 
and attitude control system with equal maximum normal accelerations. 
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